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ABSTRACT 

I study some effects of aspherical mass loss during the last stages of the asymptotic giant branch 
(AGB) on the appearance of proto planetary nebulae (PNs) and young PNs. The aspherical 
mass loss can be small scale inhomogeneities, and/or axially symmetric mass loss geometry. I first 
examine the role of the dust opacity in the optical band on the appearance of proto-PNs. I conclude 
that large optical depth will be found in proto-PNs which are descendant of AGB stars having high 
equatorial mass loss rates, which require a stellar binary companion for their existence. In these 
cases light from the central star will reach larger distances along and near the polar directions, 
leading to the appearance of an elongated reflection nebula. These proto-PNs will become bipolar- 
PNs, i.e., PNs with two lobes and an equatorial waist between them, or extreme elliptical, e.g., a 
ring, but no lobes on the two sides of the equatorial plane. I then derive the conditions for the 
enhancement of non-radial density inhomogeneity by the propagation of the ionization front at the 
early PN stages. The ionization will proceed faster in the radial direction along low density regions. 
The low density regions will be heated earlier, and they will expand due to their higher pressures, 
reducing further their densities. The opposite occurs in high density regions. The condition for this 
ionization instability to develop is that the ionization time difference between two direction at the 
same radius is longer than the sound crossing time between these two regions. This condition for 
the ionization front instability can be expressed as a condition on the mass loss rate inhomogeneity, 
i.e., its dependence on direction. 



Key words: Planetary nebulae: general — stars: AGB and post-AGB — stars: mass loss — circum- 
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1. INTRODUCTION 

High spatial resolution observations in recent years, mainly with HST, show that many proto 
planetary nebulae (PNs) and young PNs possess complex structures in their inner regions. Proto- 
PNs appear in the optical as reflection nebulae, with different optical depths along different direc- 
tions from the central star to different parts of the nebula (e.g., the Egg nebula [CRL 2688], Sahai 
et al. 1998a, b; more examples in Ueta, Meixner & Bobrowsky 2000; hereafter UMB). In most 
extreme cases, termed DUPLEX (DUst-Prominent Longitudinally-EXtended) by UMB, the light 
from the central star is completely, or almost completely, blocked in the equatorial plane. Many 
young PNs have filaments, dents, blobs and other inhomogeneities in their inner regions (e.g., Sahai 
& Trauger 1998). 

In the present paper I examine some effects of aspherical mass loss on the propagation of 
radiation during the proto-PNs and young PNs stages, which affect the structures and appearances 
of the nebulae during these stages. The aspherical mass loss can be an axisymmetrical mass loss, 
and/or local inhomogeneities in the mass loss process. Such local inhomogeneities may result from 
instabilities during the proto-PN stage (e.g., Dwarkadas Sz Balick 1998), or from the mass loss 
process itself, as expected from the model of dust formation above magnetic cool spots (Soker 2000 
and references therein). In that model, which was proposed to account for the formation of most 
elliptical PNs (Soker 2000 and references therein), the magnetic cool spots on the AGB surface are 
formed by a weak magnetic activity, and they facilitate the formation of dust. The magnetic field 
has no dynamic role. A weak dynamo activity forms more magnetic cool spots near the equatorial 
plane than closer to the poles, thus leading to axisymmetrical mass loss and the formation of 
elliptical PNs. If the dust formation occurs during the last AGB phase when mass loss rate is 
high, the dust shields the region above it from the stellar radiation (Soker 2000). This leads to 
both further dust formation in the shaded region, and, due to lower temperature and pressure, the 
convergence of the stream toward the shaded region, and the formation of a flow having a higher 
density than its surroundings. In that model the very large optical depth of the dusty envelope 
during the AGB superwind phase increases substantially the departure from spherical mass loss 
geometry, both locally and globally. This model for axisymmetrical mass loss, which is based on a 
weak dynamo activity that leads to the formation of magnetic cool spots, can explain other features 
in addition to the axisymmetrical mass loss itself (Soker & Clayton 1999). Most important, it can 
operate in very slowly rotating AGB stars, which can gain the required rotation from a planet 
companion, or else be fast rotators on the main sequence. Since the departure from spherical mass 
loss occurs only when opacity near the AGB surface is large, which means high mass loss rate, the 
model explains the higher ellipticity of the inner regions of many elliptical PNs. In extreme cases 
the inner regions are elliptical while the halos are spherical (e.g., NGC 6891, Guerrero et al. 2000). 
It can also account for local inhomogeneities, such as clumps and filaments, and it can account for 
the change in direction of the symmetry axis, which is later observed as a point symmetric elliptical 
PN. Changes in the direction of the magnetic axis is seen in the Sun (over a period of ~ 11 yrs), 
and in Earth (over a period of hundreds of thousands years) . 
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The present paper examines the role of the density inhomogeneities on radiative transfer at 
later stages: their role in the optical appearance of the reflection nebulae during the proto-PN phase 
(§2), and their role in the propagation of the ionization front in young PNs (§3). It is well established 
that the ionization front has a substantial influence on the structures of PNs (Breitschwerdt & Kahn 
1990; Mellema 1995; Mellema & Frank 1995; Chevalier 1997; Schonberner & Steffen 2000; Soker 
1998). Here (§3) I derive the condition for the density inhomogeneities to be amplified by the 
ionization front, leading to an ionization front instability. In §4 I summarize the main results. 



2. POST AGB STARS 



This section examines the optical radiation from the central star prior to ionization as it 
propagates through the dusty nebula. For the purpose of this paper it is adequate to take a simple 
form for the obscuring AGB wind. The mass loss rate per unit solid angle along a specific direction 
is m(6) = M(9)/(Air), and the wind moves radially at a constant velocity v w . In general v w may 
depend on the direction, but here I simply assume that it is uniform. I also assume that most 
of the opacity is due to the intensive mass loss episode at the termination of the AGB, i.e., the 
superwind, lasting for a time t sw . The opacity due to the low density inner regions of the post- 
AGB circumstellar matter (e.g., Schonberner & Steffen 2000) is neglected. At a time t after the 
superwind ends the density distribution is given by 

M(9)/{4irv w r 2 ) r in < r < r out 1 

very low r < r in or r > r out , J 

where the boundaries of the superwind are 

r m = v w t, and r ou t = n n + v w t sw . (2) 




The optical depth of the dusty wind is given by 

/ in K ( r ) /9 ( r ) dr 4ff Vw ( rin rout ) rin ( 1 rout ) t i 1 t + t sw ) ^ 

where the integration was performed with equation (1) for the density, and the length scale r T and 
time scale t T are defined as follows, 

= MW* ( M(9) \ / k \ ( v w \ - 1 

T ~ 4ttv w ^1O- 4 M yr" 1 J V. 150 cm 2 g~V \l0 km s' 1 J ' { ' 

and 



l ^V w =2AX 10 U-4iV yr-i J 1 150 cm 2 g-i j ^ (5) 
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The opacity of dust and gas around AGB stars is k ~ 20 cm 2 g _1 at 2fim (e.g., Jura 1986), and I 
scaled it to 150 cm 2 g _1 in the optical as it is for the ISM. At early stages, r- m <C r T , or t <C t T , the 
opacity is very large, and no light escapes at all, as is well known for many upper AGB stars. At 
very late stages r m ^> r T , or t 3> t T , and the dusty nebula is transparent. The interesting time is 
when t ~t T . 

As an example let the mass loss rate (when extended to a complete sphere) along the equatorial 
and polar directions be M(n/2) = 5 x 1O _5 M yr _1 and M(0) = 1O~ 5 M yr _1 , respectively. The 
velocity will be taken to be 10 km s _1 in all directions, and the superwind duration 1000 yrs. At 
time t = 200 yrs the optical depth of the superwind along the equatorial plane is 5, while it is 
only 1 along the polar direction. The central star is substantially attenuated along the equatorial 
plane, but only moderately so along the polar directions. Taking the equatorial flow to be slower 
will increase the effect further. 

The mass loss rate at the end of the superwind phase declines by more than two orders of 
magnitude within ~ 100 yrs (Blocker 1995; Schonberner & Steffen 2000), so any effect on a shorter 
time scale cannot be accurately treated with the assumptions made here. From equation (5) we 
learn that for the optical depth to be ^> 2, so that the central star is attenuated by a factor of ~ 10, 
for more than ~ 200 yrs the mass loss should be M{6) ^ 2x 1O~ 5 M yr _1 . This result is interesting 
since the maximum mass loss rate possible from radiation momentum transfer is M max = n s L/(cv s ), 
where L is the stellar luminosity, c is the speed of light, v w the terminal wind velocity and n s is 
the average number of times a photon is scattered within the outflowing material. For most cases 
n s £ 1 (Knapp 1986). For L = 5OOOL , v w = 10 km s _1 , and n s = 1 we find M max = 1O" 5 M yr" 1 . 
The conclusion is that the dust opacity will obscure almost completely the central star, at least 
in some directions, for systems where one or more of the following occurs: (1) luminosity is very 
high L 3> 5OOOL ; (2) the expansion velocity, v w , (probably in the equatorial plane) is very low; 
(3) mass loss rate in the equatorial plane is enhanced by dynamic effects of a binary companion. A 
very luminous central star requires a massive progenitor, which if has a non-spherical mass loss will 
probably lead to the formation of a dense equatorial flow, but not necessarily to the formation of 
a bipolar PN. In any case, a very luminous central star is required, and this is rare. Processes (2) 
and (3) are more likely to occur. We note that slowly expanding equatorial flows are found around 
several binary systems having orbital separations of ~ 1 AU (Van Winckel 1999; Van Winckel et 
al. 1998; Jura & Kahane 1999). These systems seem to form bipolar PNs (e.g., the Red Rectangle, 
Waters et al. 1998). It seems that a slowly expanding equatorial flow requires binary interaction. 
A high mass loss rate in the equatorial plane due to dynamic effects can result both from close 
companions outside the envelope (e.g., Mastrodemos & Morris 1999), or from a common envelope 
evolution as is evident from the structure of most of the 16 PNs known to have close-binary nuclei 
(orbital periods from a few hours to 16 days; Bond 2000). However, in most cases the common 
envelope evolution will not lead to the formation of a bipolar-PN but to an elliptical PN with 
high equatorial to polar density ratio, e.g., ring-like structures but without any lobes (Soker 1997). 
(Bipolar PNs are defined as axially symmetric PNs having two lobes with an 'equatorial' waist 
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between them.) 

The conclusion from this section is that high optical depth, r ^ 2, will be found in proto-PNs 
with high equatorial mass loss rates, which require their AGB progenitors to interact with stellar 
companions. The effect of this high optical depth is that light from the central star penetrates 
large distances along and near the polar directions, but not along and near the equatorial plane. 
This leads to the appearance of an elongated reflection nebula (e.g., UMB). However, not all of 
these will turn into bipolar-PNs. Some will form extreme elliptical PNs, e.g., rings without lobes 
on the two sides of the equatorial plane. UMB observe 21 reflection proto-PNs with the Hubble 
Space Telescope, and classify them into two groups. Proto-PNs with highly or completely obscured 
central stars were termed DUPLEX (DUst-Prominent Longitudinally-EXtended) reflection nebulae. 
Those with almost no obscuration are termed SOLE (Star-Obvious Low-level-Elongated) reflection 
nebula. The more elongated DUPLEX reflection nebula are mostly a result of obscuration in the 
equatorial plane. No obscuration occurs in the less elongated SOLE reflection nebula. I find that 
out of the 10 DUPLEX proto-PNs presented by UMB some do not show any signature of lobes, 
and therefore, I argue, they will later turn into extreme elliptical PNs rather than bipolar PNs. 
I further speculate that these systems may have close binary nuclei, like the 16 systems listed by 
Bond (2000). These are IRAS 19374+2359, and IRAS 23321+6545, and possibly IRAS 16342-3814, 
and IRAS 20028+3910. 



3. THE EARLY IONIZATION PHASE 



In this section I derive the condition for the enhancement of density inhomogeneities during 
the early PN phase, when the ionization front moves outward. In order to obtain an analytic 
expression, I neglect the shock preceding the ionization front (a D-type ionization front), and the 
motion of the inner boundary of the superwind inward relative to the rest of the wind, after it is 
heated by the ionization front (e.g., Breitschwerdt & Kahn 1990; Mellema 1995). I also neglect 
the compression of the inner region by the newly blown fast wind (Breitschwerdt & Kahn 1990; 
Mellema 1995; Schonberner & Steffen 2000), and any dependence on longitudes, and consider only 
axisymmetrical density profiles. I will therefore use the density profile as given in equation (1). 
However, the results are applicable to any dependence on the direction of the mass loss rate from 
the progenitor AGB star. I make some other simplifying assumptions as indicated during the 
derivation of the condition below. 

The increase in the density inhomogeneities occurs because the ionization front proceeds with 
different speeds along directions having different density profiles. Let as assume that along a 
direction 9 the density is somewhat higher that along a direction 9+ AO, because M{9) > M(0+A6), 
where the wind velocity is taken to be constant and equal along the two directions (eq. 1). Because 
of the lower density the ionization front will reach a radius r along the + AO direction earlier 
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than along the 9 direction. Let this time difference be Atj(r). The ionized region along the 9 + A9 
direction will be much hotter than the still neutral material along the 9 direction (~ 10 4 K compared 
with < 10 3 K), and its thermal pressure much higher (assuming that the density ratio is < 10), so 
it will compress the cooler region along the 9 direction. The compression proceeds in the azimuthal 
direction, from 9 + A9 to 9 with a velocity ~ c s , where c s is the adiabatic sound speed. The 
compression continues as long as the ionization front does not reach the same radius r along the 
9 direction, and it increases the density in the denser and cooler region, and decreases the density 
in the already low density region along the 9 + A9 direction. I neglect the ionization of the denser 
region along the 9 direction by the recombination radiation from the matter along the 9 + A9 
direction. The changes in the densities in the two directions will be significant only if the time for 
the compression wave to cross the distance between the two region is shorter than the time difference 
between the ionization times of the two regions Atj(r). The compression time is ~ rA9/c s . The 
condition for significant enhancement of the inhomogeneity is therefore Atj(r) > rA9/c s . Taking 
the limit of small angles and rearranging the equation, this condition reads 



Several characteristics of this instability should be noticed. First, this is not the type of 
instability considered by Breitschwerdt & Kahn (1990), who simply considered the Rayleigh- Taylor 
instability in the interface of the ionized and neutral matter. Second, if the solid angle span by 
the high density region is smaller than that of the low density region, as in a dense clump, the 
instability will result in a dense region extending radially behind the dense clump. This case was 
studied in a previous paper for a few specific cases relevant to the PN IC 4593 (Soker 1998). If it 
is the low density region that is narrower, it will be ionized first, and its density will drop further 
as it expands. This will lead to the formation of a faint region within the nebular shell. Third, 
the instability, when it starts, has a positive feedback (Soker 1998), due to the increase (decrease) 
of density of the already denser (tenuous) region, hence the ionization front will move even slower 
(faster) along this direction. Fourth, after the denser region is ionized its pressure exceeds that of 
its surroundings, and it expands and its density drops. In the present study I only examine the 
general conditions for the instability to start, and I do not follow its evolution. 

I now turn to express condition (6) in term of the density profile of the nebula. At very 
early stages the ionizing photon emission rate TV*, in photons per second, can be approximated 
by a simple linear rise with time. Using the post-AGB results presented by Blocker (1995; see 
also Schonberner & Steffen 2000) and an earlier approximation (Breitschwerdt & Kahn 1990; eq. 
1 of Soker 1998), I take the ionization to start a time t\ after the superwind ends, with a time 
dependence according to 



For the standard case I find from fig. 1 of Breitschwerdt & Kahn (1990) N = 1 x 10 47 s 1 , and 
t2 = 700 yrs, and I take the ionization to star at t\ = 1000 yrs. These numbers are very sensitive 




(7) 
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to the mass of the central star (e.g., Blocker 1995), but sufficient to illustrate, and derive the 
condition for, the enhancement of initial density inhomogeneities. For the density profile I use the 
same assumption and density profile as in the previous section (eq. 1). 

At the early stages the ionizing flux is low and the ionization front moves through a dense 
medium. At this stage most ionizing photons go to ionized recombining atoms, and only a small 
fraction of the ionizing flux reaches the ionization front and ionizes new material. I therefore 
approximate the location of the ionization front along a specific direction by equating the ionizing 
flux with the recombination rate per unit solid angle 



4tt 



/ an e nir 2 dr, (8) 

■'fin 



where a is the recombination coefficient, and n e and rn are the electron and ions number density, 
respectively. To derive the dependence of the ionization front rj(9) on time, I substitute for the 
density from equation (1) and for from equation (7), and then integrate from r Ui = v w ti to 
rf(6,ti). This gives 

rf(0,U) _ (, U-h U Y 1 ( 9 ) 



r-m V ti tF(0), 

where I defined the time scale (marked F in eq. 6 of Soker 1998) 



-i 



M\9) n e m a __ n / M{9) \ ( v w \~ N 
*F = „ 2 -77- - 700 777^771 ,„-i TTTTTZ „_i 777777^ y rs - ( 10 ) 



4vrv3 p 2 N \W- 5 M Q yr- 1 J VlOkms- 1 / \^10 47 s" 1 

More (less) massive cores have a higher (lower) mass loss rate, but the ionizing photon emission 
rate is higher (lower) as well, and ti shorter (longer). These effects may cancel, more or less, their 
mutual influence on the product t2^F m equation (9). We also note that the expression for the inner 
boundary of the superwind r- m = v w ti is not accurate, since as this region is ionized its pressure 
increases, and the hot material flows inward relative to the rest of the wind (e.g., Breitschwerdt 
& Kahn 1990; Mellema 1995). Therefore, the effective velocity in the relation r- m = v w t{ may be 
< v w . Taking rj from equation (9) for r in the instability condition (6) and taking r- m = v w ti gives 
for the instability condition 

t > !±( 1 _Uz±J*Y\ (ii) 



U c s V t 2 t F {9) 

The next step is to take the derivative of equation (9) with respect to the angle 9, at a constant 
value of rf. Rearranging equation (9), taking r- m = v w ti, and then taking the derivative gives 

d (U — ti U \ , . . 



d9 V t 2 t F {9) 

In performing the derivation we note that t\ and t% are constants, while 

dt F (6) 2 dM{9) 



d9 M(9) d9 tF ' (13) 
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so that equation (12) becomes 

/ ti - tl U ™> _t\j^_ U-h t\ \ 
f \ t 2 t F {0)M(0) t 2 t F (9) t 2 t F (6)J wi ' 1 ) 

where M' = dM/d6. Substituting rj from equation (9), and r- Ui = v w ti in equation (14) gives for t\ 

U M(9) t 2 t F \ t 2 t F ) ' 1 ' 

Elimination of t\ from equations (11) and (15) gives the instability condition on the mass loss 
inhomogeneity 

1 • - f 1 • A-) ., ' ,, • for ti >t 1 and Z>0, (16) 



M{9) c s \ t 2 t F ) 2Z(1-Zy 
where 

z(e,u)^ {ti ~ tl)u . (ir) 
t 2 t F 

Let us examine the different factors in equation (16). The sound speed of the ionized region is 
~ 12 — 14 km s _1 , depending on its temperature, while the expansion velocity is ~ 10 — 15 km s _1 . 
However, due to the slower motion of the inner boundary of the shell at early stages (e.g., Bre- 
itschwerdt & Kahn 1990; Mellema 1995), its effective value can be somewhat smaller. We therefore 
can safely take v w /c s ~ 0.5 — 1. The term [2Z(1 — Z)\~ l will reach its minimum value of 2 when 
Z = 0.5. For t\ = 1000 yrs, t 2 = 700 yrs, and tp = 700 yrs, this occurs at t, ~ 1200 yrs. By 
that time the second term on the right hand side of equation (16) is 1 + tf/t 2 tp = 3.94, and the 
instability condition becomes M'/M ^ 8. The exact minimum value of the r.h.s. of equation (16) 
(for v w = c s ) is 7.82 and it occurs at U = 1184 yrs and Z = 0.445. Taking the mass loss rate to 
be two (three) times as high as in equation (10), with all other parameters being equal, so that 
tp = 2800 yrs (6300 yrs), the minimum value of the r.h.s. of equation (16) (again, for v w = c s ) is 
4.52, (3.82) and it occurs at U = 1534 yrs and Z = 0.418 (1937 yrs and Z = 0.412). For t F > h, 
the time of maximum instability occurs at U^t\. Neglecting t\, the r.h.s. of equation (16) reads 
(v w /c s )(l + Z)/[2Z(l — Z)\. The minimum value of the r.h.s is 2.9l4(v w /c s ), and it occurs at 
Z = 2 1 / 2 - 1 = 0.414. 

The third term on the r.h.s. [2Z(1 — Z)]^ 1 will always reach a value as low as 2 when Z = 0.5. 
However, when tp <C t\ the second term (1 + tf/t 2 tp) will be very large when Z ~ 0.5, and the 
condition on the density inhomogeneity will be hard to meet. For example, when the mass loss 
rate is ~ 0.5 (~ 0.75) of that in equation (10), so that tp = 200 yrs (400 yrs), the minimum value 
of the r.h.s. of equation (16) (for v w = c s ) is 18.2 (11.0) and it occurs at t = 1062 and Z = 0.470 
(t = 1115 and Z = 0.458). 
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The main conclusion from this section is that for most elliptical PNs, density inhomogeneities 
as a result of different mass loss rates along different radial directions during the superwind phase, 
where M ^ 1O -5 M0 yr _1 , will be amplified by the propagating ionization front if M'/M ^> 4. I 
recall that the derivation of the instability condition applies to a time-independent mass loss rate, 
and hence does not apply to small blobs, but only to inhomogeneities extending to large radial 
distances. For example, for the case tF = 2800 yrs, t\ = 1000 yrs, and t2 = 700 yrs, we found above 
that maximum likelihood for the instability occurs at ti ~ 1500 yrs and Z = 0.418. From equation 
(9) we find the ionization front to be at rj = 1.7n n , hence the inhomogeneous mass loss rate lasts for 
a long time. However, the instability can start much earlier if the density inhomogeneity is larger: 
at t = 1100 yrs (1200 yrs) the condition is M'/M ;> 15 (;> 8), for which r f = 1.06r in (1.14r in ). If 
the enhanced, or reduced, mass loss rate spans an angle of 0.2 = 12°, i.e., the density inhomogeneity 
from the center of the 12° to its edge spans an angle of A9 = 0.1 = 6°, the instability condition 
is that the density will be enhanced, or reduced, by a factor of e 15Ad = 4.5 (e SAe = 2.2). We note 
that Soker (1998) finds that for a compressed tail to develop behind a dense clump, the condition 
is that the ionization front reach the clump within ~ 100 yrs of the beginning of the ionization 
(t t\ + 100 yrs), and the density enhancement be by a factor of ^ 5. The present results are 
compatible with those of Soker (1998). Dense clumps can be formed by the mass loss process 
itself, or from the interaction process of the fast and slow winds at early stages (e.g., Dwarkadas & 
Balick 1998). Dwarkadas & Balick (1998) conduct two-dimensional simulations of winds interaction, 
taking into account the evolution of the fast wind, and find that the interaction process is prone to 
instabilities, which may form clumps at early stages. Since the density enhancement is only within 
a clump, and does not extend to a large radial distance as with inhomogeneous mass loss assumed 
here, the clump density enhancement should be larger, i.e. ^ 5, to form a dense radially extended 
tail (Soker 1998). 

Images of the inner regions of several PNs observed with HST, e.g., Ml-26, He2-142, He2-138, 
and He2-131, all from Sahai & Trauger (1998), reveal many blobs, arcs and filaments with angular 
width of ~ 0.2 — 0.5 = 10 — 30°. For these inhomogeneities to be amplified by ionization, the 
density differences between different regions should have been by a factor of ^ 2 — 5 as ionization 
started. This is quite plausible in these PNs. Very narrow radially extending structures which span 
an angle of ^ 2° = 0.04, will be amplified even for mass loss rate enhancement as small as a factor 
of 1.2 - 2. 

Finally, I examine the numerical results of Mellema (1995). Mellema finds that the ionization 
front modifies the slow wind density profile along the azimuthal direction (from pole to equator) in 
his model B, but not in his model A. In his model B the density profile has a steep variation with 
the angle near the pole, where the ionization front modifies the density profile, while in model A 
the density variation with the angle is much shallower in all directions. From his density plots I find 
that within 30° from the symmetry axis (i.e., polar direction) the average value is M'/M ~ 2.5, with 
higher values nearer the symmetry axis. According to the results here, this density inhomogeneity 
can be amplified, as indeed happens in the simulation of Mellema (1995), but a shallower density 
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variation with angle, as in his model A, will not develop the ionization front instability. 



4. SUMMARY 



In the first part of the present paper (§2) I examined the role of the dust opacity in the optical 
band in the appearance of proto-PNs. This was motivated in part by the recent HST observations 
of proto-PNs by UMB, who classified the proto-PNs into two groups, SOLE and DUPLEX. The 
conclusion from this study is that a large optical depth, r ^ 2, will be found in proto-PNs with high 
equatorial mass loss rate. The high mass loss rates requires in most cases dynamic effects, probably 
from a binary companion. Such effects can be gravitational focusing by a binary companion or a 
common envelope evolution with a stellar companion. In these cases light from the central star 
will reach larger distances along and near the polar directions, leading to the appearance of an 
elongated reflection nebula. The proto-PNs will turn into bipolar-PNs, i.e., PNs with two lobes 
and an equatorial waist between them, or will become extreme elliptical, e.g., a ring, but no lobes 
on the two sides of the equatorial plane. Proto-PNs which will turn into moderate elliptical PNs, 
i.e., a small departure from sphericity, will not have high optical depth, and the light from the 
central star will not be attenuated much. Therefore, while the dust opacity near the stellar surface 
may lead to non-spherical mass loss on the upper AGB of the progenitors of moderate elliptical 
PNs (Soker 2000), it has no substantial role thereafter. 

In the second part I examined the conditions for the enhancement of non-radial density inho- 
mogeneity by the propagation of the ionization front at early stages. The ionization will proceed 
faster along low density region, which will be heated earlier than dense regions. The hot low density 
region will expand due to its higher pressure, and the density will decrease further (see Mellema 
1995 for a specific numerical simulation). The condition for this ionization instability to develop 
is that the ionization time difference between two direction at the same radius is longer than the 
sound crossing time between these two regions (eq. 6). This can be expressed as a condition on 
the mass loss variation with the direction M' = dM /d0. Assuming a constant mass loss rate with 
time, and a constant wind velocity with time and direction, this condition was derived analytically 
(eq. 16). For typical parameters of elliptical-PNs the ionization instability will increase density 
inhomogeneities when (M'/M) ^ 4. Therefore, the observed inhomogeneities in young PNs can be 
larger than the inhomogeneities of the mass loss process itself. 
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